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Quantum mechanics is a very successful and still intriguing theory, introducing 
two major counter-intuitive concepts. Wave-particle duality means that objects 
normally described as particles, such as electrons, can also behave as waves, while 
entities primarily described as waves, such as light, can also behave as particles. 
This revolutionary idea nevertheless relies on notions borrowed from classical phys¬ 
ics, either waves or particles evolving in our ordinary space-time. By contrast, en¬ 
tanglement leads to interferences between the amplitudes of multi-particle states, 
which happen in an abstract mathematical space and have no classical counter¬ 
part. This fundamental feature has been strikingly demonstrated by the violation 
of Bell’s inequalities 1 "®. There is, however, a conceptually simpler situation in 
which the interference between two-particle amplitudes entails a behaviour im¬ 
possible to describe by any classical model. It was realised in the Hong, Ou and 
Mandel (HOM) experiment®, in which two photons arriving simultaneously in the 
input channels of a beam-splitter always emerge together in one of the output 
channels. In this letter, we report on the realisation, with atoms, of a HOM ex¬ 
periment closely following the original protocol. This opens the prospect of testing 
Bell’s inequalities involving mechanical observables of massive particles, such as 
momentum, using methods inspired by quantum optics!^ 7 , with an eye on theor¬ 
ies of the quantum-to-classical transition 8 11 . Our work also demonstrates a new 
way to produce and benchmark twin-atom pair£2E5 that may be of interest for 
quantum information processing® an( j quantum simulation 1 ®. 
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A pair of entangled particles is described by a state vector that cannot be factored as a 
product of two state vectors associated with each particle. Although entanglement does not 
require that the two particles be identical®, it arises naturally in systems of indistinguishable 
particles due to the symmetrisation of the state. A remarkable illustration is the HOM ex¬ 
periment, in which two photons enter in the two input channels of a beam-splitter and one 
measures the correlation between the signals produced by photon counters placed at the two 
output channels. A joint detection at these detectors arises from two possible processes: either 
both photons are transmitted by the beam-splitter or both are reflected (Fig. [l^). If the two 
photons are indistinguishable, both processes lead to the same final quantum state and the 
probability of joint detection results from the addition of their amplitudes. Because of ele¬ 
mentary properties of the beam-splitter, these amplitudes have same modulus but opposite 


signs, thus their sum vanishes and so also the probability of joint detection (Refs. [16, 17 and 
Methods). In fact, to be fully indistinguishable, not only must the photons have the same 
energy and polarisation, but their final spatio-temporal modes must be identical. In the HOM 
experiment, it means that the two photons enter the beam-splitter in modes that are the exact 
images of each other. As a result, when measured as a function of the delay between the arrival 
times of the photons on the beam-splitter, the correlation exhibits the celebrated ‘HOM dip’, 
ideally going to zero at null delay. 


In this letter, we describe an experiment equivalent in all important respects to the HOM 
experiment, but performed with bosonic atoms instead of photons. We produce freely propagat¬ 
ing twin beams of metastable Helium 4 atom^®, which we then reflect and overlap on a beam¬ 
splitter using Bragg scattering on an optical lattice (Ref. [l9| and Fig. [TJ) . The photon counters 
after the beam-splitter are replaced by a time-resolved, multi-pixel atom-counting detectoi^, 
which enables the measurement of intensity correlations between the atom beams in well defined 
spatial and spectral regions. The temporal overlap between the atoms can be continuously 
tuned by changing the moment when the atomic beam-splitter is applied. We observe the 
HOM dip when the atoms simultaneously pass through the beam-splitter. Such a correlation 
has no explanation in terms of classical particles. In addition, a quantitative analysis of the vis¬ 
ibility of the clip also rules out any interpretation in terms of single-particle matter waves. Our 
observation must instead result from a quantum interference between multi-particle amplitudes. 
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Our experiment starts by producing a Bose-Einstein condensate (BEC) of metastable He¬ 
lium 4 atoms in the 2 3 Si,m = 1 internal state. The BEC contains 5 to 6 x 10 4 atoms and is 
confined in an elliptical optical trap with its long axis along the vertical (z) direction (Fig. [T^i) . 
The atomic cloud has radii of 58 and 5 pm along the longitudinal and transverse (_L) directions, 
respectively. A moving optical lattice, superimposed on the BEC for 300 ps, induces the scat¬ 
tering of atom pairs (hereafter referred to as twin atoms) in the longitudinal direction through 
a process analogous to spontaneous four wave mixing (Refs. (l8j 21, 221 and Methods). One 
beam, labelled a, has a free-space velocity v z ~ 12.1 cms' 1 in the laboratory frame of reference 
and the other beam, labelled b , has a velocity v z ~ 7.0 cm s^ 1 (Fig. [lj3,c) . The twin atom 
beams clearly appear in the velocity distribution of the atoms, which is displayed in Fig. [2] 
The visible difference in population between the beams probably results from secondary scat¬ 
tering processes in the optical lattice, leading to the decay over time of the quasi-momentum 
stated 8 . After the optical lattice has been switched off (time t±), the twin atoms propagate 
in the optical trap for 200 ps. At this moment, the trap itself is switched off and the atoms 
are transferred to the magnetically insensitive rn — 0 internal state by a two-photon Raman 
transition (Methods). 

From here on, the atoms evolve under the influence of gravity and continue to move apart 
(Fig. [I)o). At time t 2 = f i + 500 ps, we deflect the beams using Bragg diffraction on a second 
optical lattice, so as to make them converge. In the centre-of-mass frame of reference, this 
deflection reduces to a simple specular reflection (Fig. [lj: and Methods). At time f 3 ~ 2 t 2 — t\, 
we apply the same diffraction lattice for half the amount of time in order to realise a beam- 
splitting operation on the crossing atom beams. Changing the time t 3 allows us to tune the 
degree of temporal overlap between the twin atoms. Fig. 03 shows the atomic trajectories 
in the centre-of-mass frame of reference and reveals the close analogy with a photonic HOM 
experiment. The atoms end their fall on a micro-channel plate detector located 45 cm below 
the position of the initial BEC and we record the time and transverse position of each atomic 
impact with a detection efficiency p ~ 25% (Methods). The time of flight to the detector is 
approximately 300 ms, long enough that the recorded signal yields the three components of the 
atomic velocity. By collecting data from several hundred repetitions of the experiment under 
the same conditions, we are able to reconstruct all desired atom number correlations within 
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variable integration volumes of extent Av z x Av\. These volumes play a similar role to that 
of the spatial and spectral filters in the HOM experiment and can be adjusted to erase the 
information that could allow tracing back the origin of a detected particle to one of the input 
channels. 

The HOM dip should appear in the cross-correlation between the detection signals in the 
output channels of the beam-splitter (Ref. [l6| and Methods): 
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Here, a v and fit denote the annihilation and creation operators of an atom with three-dimensional 
velocity v, respectively, (•) stands for the quantum and statistical average and V c ,d designate 
the integration volumes centred on the output atom beams c and d (Fig. [l^). We have measured 
this correlation as a function of the duration of propagation r = t 3 — f 2 between the mirror and 
the beam-splitter (Fig. [3]) and for various integration volumes (see supplementary material). 
We observe a marked reduction of the correlation when r is equal to the duration of propaga¬ 
tion from the source to the mirror (£3 — £2 — £2 — £ 1 ) and for small enough integration volumes, 
corresponding to a full overlap of the atomic wave-packets on the beam-splitter. Fitting the 
data with an empirical Gaussian profile yields a visibility: 

v= max T (r) - min T G { c 2 J(t 
max T G ( c 2 J(t) 

where the number in parenthesis stands for the 68 % confidence interval. As we shrink the integ¬ 
ration volumes, we observe that the dip visibility first increases and then reaches a saturation 
value, as is expected when the integration volumes become smaller than the elementary atomic 
modes. The data displayed in Fig. [ 3 ] were obtained for Av z = 0.3cms _1 and Av ± = 0.5 cm s ', 
which maximises the reduction of the correlation while preserving a statistically significant 
number of detection events (see supplementary material). 

The dip in the cross-correlation function cannot be explained in terms of classical particles, 
for which we would have no correlation at all between the detections in the output channels. 
When the atoms are viewed as waves however, demonstrating the quantum origin of the dip 
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necessitates a deeper analysis. The reason is that two waves can interfere at a beam-splitter 
and give rise to an intensity imbalance between the output ports. If, in addition, the coherence 
time of the waves is finite, the cross-correlation can display a dip similar to the one observed in 
our experiment. But once averaged over the phase difference between the beams, the visibility 
is bounded from above and cannot exceed 0.5 (Refs. 123, 24] and Methods). In our experiment, 
this phase difference is randomised by the shot-to-shot fluctuations of the relative phase between 
the laser beams used for Bragg diffraction (Methods). Since our measured visibility exceeds 
the limit for waves by two standard deviations, we can safely rule out any interpretation of our 
observation in terms of interference between two ‘classical’ matter waves or, in other words, 
between two ordinary wave functions describing each of the two particles separately. 

Two contributions may be responsible for the non-zero value of the correlation function at 
the centre of the dip: the detected particles may not be fully indistinguishable and the number 
of particles contained in the integration volume may exceed unity for each beam (Refs. 123, 
and Methods). The effect of the atom number distribution can be quantified by measuring the 
intensity correlations of the twin atom beams upstream of the beam-splitter (Fig. [IJ 3 ), which 
bound the visibility of the dip through the relation: 
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where G„ 2 J , and G ( ^j are defined according to Eq. [l] (Ref. |24j and Methods). Here, one 
immediately sees that the finite probability of having more than one atom in the input channels 
will lead to finite values of the auto-correlations Gad , G hb and therefore to a reduced visibility. 
We have performed the measurement of these correlations following the same experimental 
procedure as before, except that we did not apply the mirror and beam-splitter. We find non¬ 
zero values G^aa = 0.016(5), G^f = 0.047(9), and G^j = 0.048(7), yielding V max = 0.60(10), 
where the uncertainty is the standard deviation of the statistical ensemble. Because of the good 
agreement with the measured value, we conclude that the atom number distribution in the input 
channels entirely accounts for the visibility of the HOM dip. For the present experiment we 
estimate the average number of incident atoms to be 0.5(1) in V a and 0.8(2) in V&, corresponding 
to a ratio of the probability for having two atoms to that for having one atom of 0.25(5) and 
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0.40(10), respectively (Methods). Achieving much smaller values is possible, for instance by 
reducing the pair production rate, but at the cost of much lower counting statistics. 

Although multi-particle interferences can be observed with particles emitted or prepared 
independently 13 26 29 , twin particle sources are at the heart of many protocols for quantum 
information processing 14 and quantum simulation 15 . The good visibility of the HOM dip in 
our experiment demonstrates that our twin atom source produces beams which have highly 
correlated populations and are well mode matched. This is an important achievement in itself, 
which may have the same impact for quantum atom optics as the development of twin photon 
sources using non-linear crystals had for quantum optics (see for instance Ref. |30|). 
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Figure 1 | Schematic of the experiment, a, A Bose-Einstein condensate (BEC, 
white oval) of metastable Helium 4 atoms is trapped in an elongated optical trap (red 
shaded area). A moving optical lattice, here depicted in blue, is superimposed on the 
BEC and triggers the scattering of atom pairs along the z-axis. When this lattice and the 
trap are switched off, the atoms fall towards a micro-channel plate detector located 45 cm 
below the initial position of the BEC (drawing not to scale), b, The time diagram shows 
the evolution of the twin atom vertical coordinates (blue and red lines). Between t\ and 
t 2 , t 2 and t 3 , and after t 3 , the atoms move under the sole influence of gravity (drawing 
not to scale). At t 2 , the twin atom velocities are swapped using Bragg diffraction on an 
optical lattice. At time f 3 , when the atomic trajectories cross again, the same lattice is 
applied for half the amount of time in order to realise a beam-splitter, c, In the centre- 
of-mass frame of reference, the trajectories of the atoms resemble those of the photons 
in the Hong-Ou-Mandel experiment. A joint detection arises either when both atoms 
are transmitted through the beam-slitter (upper panel) or when both are reflected (lower 
panel). If the two particles are indistinguishable, these processes end in the same final 
quantum state and the probability of joint detection results from the addition of their 
amplitudes. For bosons these amplitudes have same modulus but opposite signs, thus 
their sum vanishes and so also the probability of joint detection. 
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Figure 2 | Velocity distribution of the twin atoms, a, Two-dimensional velocity 
distribution of the twin atom beams emitted by the source. The red shaded area, drawn 
here only for the lower beam, labelled b in Fig. [TJ) and c, shows the integration volume used 
for computing the correlation function displayed in Fig. [3} The distribution corresponds to 
an average over about 1,100 measurements and is not corrected for the limited detection 
efficiency. The velocities are given relative to the laboratory frame of reference. The 
size of each pixel is 0.25cms _1 in the transverse directions (x and y) and 0.15cms _1 in 
the longitudinal ( z ) direction and an integration over 2 pixels is performed along the y 
direction. b,c, Cuts of the two-dimensional velocity distribution through the centre of 
the lower beam along the longitudinal (b) and transverse (c) directions. The data points 
result from the average over 2 pixels along the direction perpendicular to the cut. The 
full width at half-maximum of the distribution, obtained from a Gaussian fit, is about 
1.4cms _1 along both the longitudinal and transverse directions. The red shaded area 
again shows the integration volume. 
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Figure 3 | HOM dip in the cross-correlation function. The correlation G ( f ( j between 
the output ports of the beam-splitter, defined in Eq. [lj was measured as a function of 
the duration of propagation r = £3 — £2 between the mirror and the beam splitter. The 
HOM dip is directly visible as a marked reduction of the correlation when r equals the 
duration of propagation t 2 — ti ~ 500 ps between the source and the mirror, corresponding 
to symmetric paths between the source and the beam-splitter, i.e. when one cannot 
distinguish between the two diagrams of Fig. [lj:. A Gaussian fit (blue line) precisely 
locates the dip at r = 550(50) ps, with a full-width at half-maximum of 150(40) ps, where 
the uncertainty corresponds to the 68 % confidence interval. The measured visibility 
is V = 0.65(7). It is two standard deviations beyond the classical-to-quantum threshold 
represented by the red shaded area, which takes into account the experimental uncertainty 
over the background correlation value. Each data point was obtained from an average 
over about 500 to 1,400 repetitions of the experiment. Error bars denote the standard 
deviation of the statistical ensemble. The mean detected atom number was constant over 
the range of values of r displayed here (see supplementary material). 
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Methods 


Twin atom source 

The twin atom beams result from a scattering process between pairs of atoms from the BEC 
occurring when the gas is placed in a moving one-dimensional optical lattice. The experimental 
set-up has been described in Ref. [18]. The lattice is formed by two laser beams derived from 
the same source emitting at the wavelength A = 1,064 nm. In contrast to our previous work, 
the axis of the optical lattice was now precisely aligned with the long axis of the optical trap 
confining the atoms. The laser beams intersect with an angle of 9 — 166°, their frequency 
difference is set to 100.5 kHz and the lattice depth to 0.4f? rec . This constrains the longitudinal 
wave-vector of the twin atoms to the values k Zt<1 = 0.75 kr ec and k z j> = 1.30 £^ ec in order to 
fulfil the conservation of quasi-momentum and energy in the frame co-propagating with the 
lattice. Here, k vec = 27rsin(#/2)/A is the recoil wave-vector along the longitudinal axis gained 
upon absorption of a photon from a lattice laser and E rec = h 2 k 2 ec /2m is the associated kinetic 
energy, with h the reduced Planck constant and m = 6.64 x 10“ 27 kg the mass of an Helium 4 
atom. The observed velocities of the twin atom beams coincide with the expected values above, 
using the relation v = hk/m. The optical lattice is turned on and off adiabatically so as to 
avoid diffraction of the atoms during this phase of the experiment. 

Transfer to the magnetically insensitive internal state 

Transfer to the m — 0 state after the optical trap has been switched off is made necessary 
by the presence of stray magnetic fields in the vacuum chamber that otherwise would lead 
to a severe deformation of the atomic distribution during the long free fall. The transfer is 
achieved by introducing a two-photon coupling between the m — 1 state, in which the atoms 
are initially, and the m — 0 state using two laser beams derived from a single source emitting at 
1,083 nm and detuned by 600 MHz from the 2 3 Si to 2 3 P 0 transition. The frequency difference 
of the laser beams is chirped across the two-photon resonance so as to realise an adiabatic fast 
passage transition (the frequency change is 300 kHz in 300 ps). We have measured the fraction 
of transferred atoms to be 94%. The remaining 6% stay in the m — 1 state and are pushed 
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away from the integration volumes by the stray magnetic field gradients. 

Atomic mirror and beam-splitter 

The mirror and beam-splitter are both implemented using Bragg scattering on a second optical 
lattice. This effect can be seen as a momentum exchange between the atoms and the laser beams 
forming the lattice, a photon being coherently absorbed from one beam and emitted into the 
other. In our experiment, the laser beams forming the lattice have a waist of 1 mm and are 
detuned by 600 MHz from the 2 3 Si to 2 3 P 0 transition (they are derived from the same source as 
the beams used for the Raman transfer). In order to fulfil the Bragg resonance condition for the 
atom beams, the laser beams are made to intersect at an angle of 32° and the frequency of one 
of the beams is shifted by 53.4 kHz. In addition to this fixed frequency difference, a frequency 
chirp is performed to compensate for the acceleration of the atoms during their free fall. The 
interaction time between the atoms and the optical lattice was 100 ps for the mirror operation 
(7r-pulse) and 50 ps for the beam-splitter operation (7r/2-pulse). The resonance condition for 
the momentum state transfer is satisfied by all atoms in the twin beams but only pairs of 
states with a well defined momentum difference are coupled with each other. We measured the 
reflectivity of the mirror and the transmittance of the beam-splitter to be 0.95(2) and 0.49(2), 
respectively. Spontaneous scattering of photons by the atoms was negligible. 


Detection efficiency 


Our experiment relies on the ability to detect the atoms individually. The detection efficiency 
is an essential parameter for achieving good signal to noise ratios, although it does not directly 
influence the visibility of the HOM dip. Our most recent estimate of the detection efficiency 
relies on the measurement of the variance of the atom number difference between the twin 


beams. For this we use the same procedure as described in Ref. 18 , but with an integration 


volume that includes the entire velocity distribution of each beam. We ford a normalised 
variance of 0.75(5), well below the Poissonian floor. Since for perfectly correlated twin beams 
the measured variance would be 1 — rj, we attribute the lower limit of 25(5) % to our detection 
efficiency. This value for p is a factor of about 2 larger than the lower bound quoted in Ref. 
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The difference is due to the change of method employed for transferring the atoms from the 
m — 1 to the m — 0 state after the optical trap has been switched off. We previously used 
a radio-frequency transfer with roughly 50 % efficiency whereas the current optical Raman 
transfer has close to 100 % efficiency. 

Distribution of the number of incident atoms 

We have estimated the average number of incident atoms in each input channel of the beams- 
splitter, n a and rib, by analysing the distribution of detected atoms in the integration volumes 
V a and Vb- We fitted these distributions by assuming an empirical Poissonian law for the 
distribution of incident atoms and taking into account the independently calibrated detection 
efficiency. The values of n a and n& given in the main text are the mean values of the Poissonian 
distributions that best fit the data. The probabilities for having one or two atoms in each of 
the input channels of the beam splitter was obtained from the same analysis. The uncertainty 
on these numbers mostly stems from the uncertainty on the detection efficiency. 


The HOM effect 


The HOM effect appears in the correlator (a^ af rf a Vrf a Vc ) of Eq. [lj The simplest way to calculate 
such a correlator is to transform the operators and the state vector back in the input space 
before the beam-splitter and to use the Heisenberg picture. The transformation matrix between 
the operators a Vc (t 3 ), a Vd (t 3 ) and a Va (t 3 ), a Vb (t 3 ) can be worked out from first principles. For 
the Bragg beam-splitter, and using a Rabi two-state formalism, we find: 
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where 0 is the relative phase between the laser beams forming the optical lattice. In the ideal 
case of an input state with exactly one atom in each channel, |l Va , l Vi) ), we therefore obtain: 
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meaning that the probability of joint detection is strictly zero. The detailed calculation above 
makes clear that the perfect destructive interference between the two-particle state amplitudes 
associated with the two diagrams of Fig. [lj^ is at the heart of the HOM effect. By contrast, input 
states containing more than one atom per channel are transformed into a sum of orthogonal 
states and the interference can only be partial. Taking |2 Vo ,2 v& ), for instance, yields: 
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Finally, we note that losses in one of the incident beams, for instance beam a, can be mod¬ 
elled by a fictitious beam-splitter with a transmission coefficient T. In the above calculation, 
these losses would therefore only manifest by an additional factor y/T in front of every operator 
a Va , leaving unaffected the destructive interference that gives rise to the HOM effect. 

Visibility of the HOM dip 

A slightly less general form of Eq. [3] has been derived in Ref. |24| assuming a two-mode squeezed 
state as an input state. The same calculation can be performed for an arbitrary input state. 
Leaving aside the integration over the velocity distribution, we find that the cross-correlation 
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for indistinguishable particles can be expressed as: 
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whereas that of distinguishable particles reads 

G S’L. = j( G S + G » +2G2>) ■ (13) 

Here, the correlators appearing in the right-hand side are taken at time t±, that is immediately 
after the the atom beams have been produced. The term A corresponds to an interference 
between single-particle matter waves, ft depends on both the relative phase between the atom 
beams and the relative phase between the laser beams used for Bragg diffraction. The latter 
is counted once for the atomic mirror and once for the atomic beam-splitter. Twin beams 
with perfect correlations in their population would have a fully random relative phase. In our 
experiment however, the population imbalance between the atom beams could entail a residual 
phase coherence. Instead, the relative phase between the laser beams was left uncontrolled and 
its value was randomly distributed between two repetitions of the experiment. As a result, 
the term A must average to zero and the visibility of the HOM dip be given by Eq. [3j as 
observed in the experiment. Following Ref. |24|, we also note that Eq. [3] entails the ultimate 
bound for waves interfering on the beam-splitter: because waves must fulhl the Cauchy-Schwarz 
inequality, < \J , the visibility of the classical clip cannot exceed 0.5. 

The above results holds true for a finite integration over the atomic velocity distribution 
if two conditions are met: (i) It must remain impossible to distinguish the atoms entering 
the beam-splitter through channel a from the atoms entering through channel b once they have 
exited the beam-splitter; (ii) The transformation matrix of the beam-splitter must keep the same 
form after integration. In our experiment, the second condition is naturally satisfied because 
the Bragg diffraction only couples atomic states with a well defined momentum difference and 
we fulhl the first condition by reducing the integration volume as much as it is necessary. 
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Supplementary material 


Optimization of the Hong-Ou-Mandel dip 

The visibility of the Hong-Ou-Mandel dip is plotted in Fig. [ST] as a function of the longitudinal 
(a) and transverse (b) integration volume. The red dots identify the integration volume used 
in Fig. [3] of the main text and correspond to a compromise between signal to noise ratio and 
visibility amplitude. As we shrink the integration volumes, the dip visibility first increases and 
then reaches a saturation value, meaning that the integration volume becomes smaller than 
the elementary atomic moded 32-34 . Reducing further the integration volume only leads to an 
increase of the statistical uncertainty. 

The visibility V is obtained by fitting the cross-correlation function (r) measured in the 
experiment with the empirical function: 

f( T ) = G b 2 g (l - Rexp ( - (t - t 0 ) 2 /2 (x 2 )) , 

( 2 ) 

where the background correlation G ^ , the center of the dip r 0 and the width of the dip a are 
all left as free parameters. 
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Figure SI | Visibility as a function of the integration volume, a, Visibility as a 
function of the longitudinal integration interval Av~. The transverse integration interval 
is kept constant at Ar?j_ = 0.48cms _1 . b, Visibility as a function of the transverse 
integration interval Ar?_|_. The longitudinal integration interval is kept constant at Av z = 
0.28 cm s -1 . The red points mark the values discussed in the main text. 


Stability of the atom number in the output ports 


The mean detected atom number in the output ports c and d is plotted as function of r in 
Fig. |S2j The mean atom number is constant as function of r within the statistical uncertainty. 
To easily compare the atom number fluctuations to the variation of the cross-correlation across 
the HOM dip, the product of the averaged populations {n c )-{nd) and the cross-correlation 
are displayed together as a function of r in Fig. [S2f . In contrast to the cross-correlation, it is 
impossible to identify a marked reduction of (n c )-(rid) around r = 550 ps. This confirms our 
interpretation of the dip as a destructive two-particle interference. 
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Figure S2 | Stability of the output population over propagation duration, a, 

Averaged atom number detected in V c as a function of the propagation time r. The mean 
value of n c is 0.20 with a standard deviation of 0.01. b, Averaged atom number detected 
in Vd as a function of the propagation time r. The mean value of rid is 0.19 with a 
standard deviation of 0.01. c, The cross-correlation between the output ports c and d 
(solid blue circles), corresponding to the HOM dip, is compared to the product of the 
average density populations (n c ) • (rid) (open gray circles). The product of the averaged 
population is constant while the cross correlation exhibits a dip around r = 550 ps. 
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